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1.   INTRODUCTION 

The  model  presented  in  this  paper  is  an  event  store  computer 
simulation  of  the  interactions  between  surface-to-air  missile 
systems  and  aircraft  in  a  non-jamming  environment  and  over  flat 
terrain.  The  purpose  of  the  model  is  to  test  the  general  dis- 
position of  missile  areas  and  the  associated  missile  system  reaction 
times  against  an  aircraft  attack.  The  model  is  a  probabilistic 
simulation  in  that  the  success  or  failure  of  any  probabilistic 
event  is  determined  in  the  model  by  comparing  the  numerical  value 
assigned  to  the  probability  of  success  or  failure  to  a  program 
generated  random  number.   It  is  the  case  in  this  model  that  all 
numbers  so  used  are  uniformly  distributed.  A  complete  description 
of  the  model  is  included  in  the  following  sections. 

V7  2.  PLAYING  AREA 

The  playing  area  for  the  model  is  that  portion  of  a  circle 
whose  center,  radius  and  central  angle  are  inputs.  The  structure 
of  the  computer  program  is  such  that  the  playing  area  must  be 
located  entirely  within  the  first  quadrant  of  the  X,Y  plane.  The 
numerical  restrictions  within  the  computer  program  are  such  that 
the  central  angle  and  radius  must  be  less  than  180  degrees  and 
1000  miles  respectively.  The  playing  area  must  also  be  oriented 
such  that  the  bisector  of  the  central  angle  is  parallel  to  the 
X-axis  while  the  central  angle  opens  in  the  positive  X-direction. 
The  playing  area  is  illustrated  in  Figure  1. 


*>x,gy) 


(0,0) 


CGX,G/) 
~RCC 
9 


Coor<//*a^es  of  Center 
Hdius  of  7^/cLYna  Area 
Centra/     Ay<*/& 


A/oie ; 


0    <   GX  <    GX  +  7(CC- Cos  d/2. 

o  <  Gy  -   7?cc  •  Sin  9/l   <    6y 


7?CC 

e 


/OQO    ~*nJe,S 

/So      Je&rees 


ftquve    / 


V    \ 


f 


y 


3.  OFFENSE 

The  offense  consists  of  as  many  as  twenty   aircraft. 
These  aircraft  fly  through  the  playing  area  in  an  attempt  to 
penetrate  the  missile  defenses.  The  entry  points  into  the  playing 
area  for  the  aircraft  are  generated  by  the  computer  program  by 
assuming  the  entry  points  are  uniformly  distributed  along  the  arc 
of  the  circle  defined  by  the  radius  of  the  playing  area  and  the 
Central  angle.  The  flight  path  for  each  aircraft  after  it  enters 
the  playing  area  is  to  fly  straight  toward  the  center,  (GX,GY)  . 

The  spacing  time  between  aircraft  and  the  speeds  and  altitudes  ^ 

of  aircraft  are  also  generated  by  the  computer  program.  These 

/ 
values  are  assumed  by  the  model  to  be  uniformly  distributed  between 

their  respective  minimum  and  maximum  values.  These  minimum  and      * 

maximum  values  are  inputs  to  the  model. 

The  aircraft  in  the  model  play  a  passive  role  and  serve  only    / 

as  the  set  of  stimuli  needed  to  cause  the  missile  systems  to  act. 

These  aircraft  do  not  defend  themselves  against  missile  attack  nor 

do  they  attack  the  missile  areas.  Under  the  assumption  that  the 

center  of  the  playing  area  represents  the  bomb  release  line  of  the 

aircraft  for  their  respective  targets,  an  offensive  wave  of  aircraft 

is  considered  successful  if  at  least  one  aircraft  reaches  the  center.  J 

4.  DEFENSE 

The  defense  consists  of  as  many  as  three  missile  areas  with 
their  associated  missile  systems.  These  missile  areas  may  be 
located  anywhere  in  the  first  quadrant.  They  need  not  be  located 
within  the  playing  area;  however,  since  only  the  results  of 


V 


r^AKM 


^4uAJ 


interactions  occurring  within  the  playing  area  are  considered  in 
the  model* the  sphere  of  influence  of  the  missile  area  must  include 
some  portion  of  the  playing  area  in  order  for  the  missile  areas 
to  exert  any  effect  on  the  simulation  results. 

Associated  with  each  missile  area  are  the  parameters  needed 
to  describe  its  missile  system.  The  values  of  these  parameters 
are  inputs  to  the  model  and  the  parameters  are: 

(1)  Search  radar  maximum  range 

(2)  Missile  maximum  range 

(3)  Missile  average  speed 

(4)  The  number  of  tracking  radars 

(5)  The  number  of  missile  launchers 

(6)  Maximum  and  minimum  time  required  to  reload  a  launcher 

(7)  Maximum  and  minimum  time  required  to  acquire  a  target 

on  the  tracking  radar 

(8)  Maximum  and  minimum  time  required  to  assess  a  target 

after  missile  intercept 

(9)  Missile  single-salvo  kill  probability 

It  is  an  assumption  of  the  model  that  all  aircraft  are  observed 
by  all  missile  areas  subject  to  the  aircraft  radar  horizon  and  the 
missile  area  search  radar  maximum  range.  It  is  also  the  case  that 
in  order  to  fire  a  missile,  or.  salvo,  at  an  aircraft: 

(1)  The  aircraft  must  be  observed  at  the  time  of  fire. 

(2)  A  missile  launcher  must  be  loaded. 


(3)  A  tracking  radar  must  be  free  in  order  to  be  used  for 
full  course  missile  guidance. 

(4)  The  intercept  point  must  be  within  the  missile  maximum 
range  circle. 

(5)  The  aircraft  must  not  be  past  the  point  of  closest 
approach  to  the  missile  area  at  the  time  of  fire. 

The  model  does  not  include  altitude  or  minimum  range  restrictions 
bn  the  missile. 

The  significant  time  delays  inherent  to  the  missile  systems 
and  included  in  the  model  are  seen  to  be: 

(1)  Reload  time:  The  amount  of  time  required  to  reload  a 
missile  launcher. 

(2)  Acquisition  time:  The  amount  of  time  required,  once  an 
aircraft  is  observed  on  the  search  radar,  to  transfer  the  aircraft 
as  a  target  to  an  available  tracking  radar. 

(3)  Assessment  time:  The  amount  of  time  the  tracking  radar 
must  remain  trained  on  the  target  after  missile  intercept  in  order 
for  the  results  of  the  intercept  to  be  observed. 

All  of  these  times  when  used  in  the  model  are  assumed  to  be 
uniformly  distributed  between  their  maximum  and  minimum  values. 

The  firing  doctrine  for  a  missile  system  is  shoot-look-shoot 
at  all  available  aircraft.  That  is,  when  a  missile  area  has 
launched  a  salvo  against  a  target  no  new  salvos  against  that  target 
will  be  launched  from  that  missile  area  until  that  salvo  has 
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intercepted  the  target  and  the  results  of  the  intercept  have  been 
assessed.  The  aircraft  are  selected  as  targets,  within  the 
missile  launcher  and  tracking  radar  numerical  restrictions,  on  a 
first-come  first-served  basis. 

Many  of  the  built  in  characteristics  of  the  missile  system 
can  be  overridden  by  a  proper  selection  of  the  input  parameter 
values.  For  instance,  by  selecting  a  very  high  value,  say  99, 
for  the  number  of  tracking  radars  at  a  missile  area  the  availability 
of  tracking  radars  will  not  effect  the  simulation  results,  i.e., 
the  missile  would  be  considered  as  a  non-beam  rider.  This  example 
is  based  upon  the  firing  doctrine  of  the  missile  area  and  the 
number  of  aircraft  in  the  model. 

An  illustration  of  the  playing  area  with  a  typical  missile 
area  and  aircraft  flight  path  is  included  as  Figure  2* 

5.  GAME  DOCTRINE 

With  the  input  parameter  values  assigned  the  model  considers 
the  interactions  that  occur  in  the  playing  area  between  the  missile 
systems  and  aircraft.  For  the  given  set  of  defensive  and 
offensive  parameters  the  required  number  of  aircraft  will  enter 
the  playing  area  at  points,  times,  speeds  and  altitudes  generated 
by  the  computer  program.  This  set  of  aircraft  will  then  proceed 
directly  toward  the  center,  (GX,GY) ,  passing  through  the  missile 
defenses. 
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One  complete  pass  through  the  computer  simulation  with  one 
set  of  aircraft  is  referred  to  as  a  replication.  To  generate  data 
for  statistical  purposes,  at  the  completion  of  a  replication  the 
computer  program  will  generate  a  new  set  of  aircraft  and  using  the 
same  set  of  input  values  will  produce  another  replication.  The 
desired  number  of  replications  is  an  input  value  and  must  be  less 
than  twenty-one.  An  entire  set  of  replications  for  a  given  number 
of  aircraft  is  referred  to  as  a  run.  The  model  output  then  consists 
of  any  of  the  following  forms  of  output: 

(1)  Battle  History:  An  event  history  of  each  replication 
containing  the  generated  events  of  the  battle  in  the  order  in  which 
the  events  occur  and  are  generated. 

(2)  Standard:  A  compilation  of  each  replication  containing  all 
aircraft  initial  conditions  and  the  number  of  salvos  fired  by  each 
missile  area  at  each  aircraft  and  the  identification  of  the  missile 
area  responsible  for  killing  each  aircraft. 

(3)  Summary:  A  summary  of  information,  by  totals  with  respect 
to  replication,  for  each  run  including  the  sample  mean,  variance 
and  standard  deviation  of  all  totals  presented. 

The  computer  program  will  then  make  as  many  runs  as  desired 
by  the  model  user  with  an  increased  number  of  aircraft  for  each  run. 
The  number  of  aircraft  in  the  first  run,  the  increment,  for  the 
number  of  aircraft  in  each  new  run,  and  the  number  of  runs  are 
input  values.  Of  interest  is  that  each  new  run  is  considered  by 
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the  model  to  be  an  extension  of  the  previous  run,  that  is,  if  run 

V 

three  contained  seven  aircraft  and  run  four  is  to  contain  nine 
aircraft,  then  for  all  replications  in  run  four  the  first  seven 
aircraft  will  have  entry  points,  altitudes,  speeds  and  times 
identical  to  those  replications  in  run  three,  etc.  It  is  also 
the  case  that  the  random  numbers  used  in  the  replications  of  a 
run  in  order  to  determine  the  outcome  of  probabilistic  events  will 
be  used  again  in  the  replications  of  a  new  run,  that  is,  the 
starting  point  in  the  sequence  of  random  numbers  used  for  repli- 
cation two  of  run  three  will  be  the  same  as  the  starting  point  in 
the  sequence  for  replication  two  of  run  four.   In  this  manner  it 
is  hoped  that  any  changes  in  the  results  between  runs  three  and 
four  can  be  attributed  to  the  increase  in  the  number  of  aircraft 
rather  than  to  the  deviations  of  the  sets  of  random  numbers  used. 
The  sequences  of  random  numbers  used  can  also  be  controlled 
by  the  model  user  in  order  to  reduce  to  a  minimum  the  effect  of 
gross  random  number  fluctuations  on  the  simulation  results.  The 
controlling  input  values  are  the  initial  starting  point  in  the 
random  number  sequence  for  replication  one  of  run  one  and  an 
increment  used  to  determine  the  starting  points  of  all  succeeding 
replications.  These  values  can  be  selected  such  that  the  random 
number  sequence  used  in  a  replication  will  be  generated  so  that 
it  overlaps  the  sequence  used  in  the  previous  replication.  This 
is  illustrated  as  follows.   If 
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represents   the  sequence  of  random  numbers  used   in  replication 
one  of  run  one  where  RNl,l  and  RNl,L  are  the  first  and   last 
numbers  used  then  the  sequence  for  replication  two  can  be  made  to 
overlap  this  sequence  by  choosing  an  increment,  DRN,    such  that 
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This  overlap  will  then  continue   throughout  the  replications  as 
illustrated: 
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The  numerical  value  assigned  to  DRN  is  not  the  particular  value 
of  any  random  number  but  the  number  of  random  numbers  between 
RNl,l  and  RN2,1. 
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When  a  new  run  is  made  the  computer  program  random  number 
generator  is  reset,  as  indicated  in  the  previous  paragraph,  to 
begin  the  sequence  for  replication  one  of  the  new  run  at  RNl,l. 
The  same  value  of  DRN  is  again  used  and  in  this  way  the  sequences 
of  random  numbers  used  for  the  new  replications  will  correspond 
to  the  sequences  used  for  the  replications  of  the  previous  run. 

The  computer  program  is  constructed  such  that  either  or  both 
of  two  missile  firing  procedures  may  be  used.  These  procedures 
are  referred  to  as  uncoordinated  and  coordinated  and  the  procedure 
used  in  the  simulation  is  determined  by  the  model  user  as  an  input 
to  the  model.  The  uncoordinated  missile  firing  procedure  allows 
all  missile  areas  in  the  simulation  to  fire  missiles  at  all  aircraft 
that  can  possibly  be  fired  upon  while  the  coordinated  missile 
firing  procedure  allows  a  missile  area  to  fire  missiles  at  an  aircraft 
only  if  no  other  missile  area  is  currently  engaging  that  aircraft. 
For  any  given  set  of  input  parameters  the  computer  program  will 
make  the  necessary  replications  and  runs  of  the  simulation  using 
either  or  both  of  these  procedures.  Of  interest  is  that  when  the 
model  user  elects  to  employ  both  procedures,  they  are  not  intermixed 
in  the  simulation  but  are  run  separately  and  simulation  results  are 
presented  so  that  comparisons  can  be  made  as  a  function  of  the 
procedure  used.   It  is  also  the  case  that  when  both  missile  firing 
procedures  are  used  the  same  sets  of  aircraft  and  sequences  of 
random  numbers  are  used  in  the  corresponding  replications  and  runs 
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of  the  simulation.  This  is  illustrated  in  Figure  3  which  displays 
the  overall  simulation  structure. 

6.  EVENTS 

As  mentioned  earlier  the  model  is  an  event  store  computer 
simulation,  i.e.,  all  actions  that  are  to  occur  in  the  simulation 
are  dynamically  generated  by  the  computer  program  as  a  result  of 
previous  simulation  actions  and  are  listed  chronologically  in  an 
Event  Store  List.  Each  of  the  actions  included  in  the  simulation 
assumes  the  form  of  a  computer  program  subroutine,  called  an  event, 
and  the  information  pertaining  to  the  action  on  the  Event  Store  List 
is  the  information  needed  to  execute  the  proper  subroutine.  There 
are  only  four  major  actions  included  in  the  simulation  as  events 
and  these  events  are: 

(1)  Fire  Missile  Salvo 

(2)  Missile  Intercept 

(3)  Reload  Missile  Launcher 

(4)  Free  the  Tracking  Radar  from  an  Intercepted  Target 

Each  of  the  computer  program  subroutines  representing  these  events 
uses  as  input  parameters  the  following  information: 

(1)  Time  event  is  to  occur 

(2)  Identification  of  Event 

(3)  Identification  of  Aircraft 

(4)  Identification  of  Missile  Area 

A  fifth  event,  Change  Engagement  Status,  is  also  part  of  the 
simulation  but  does  not  represent  an  interaction  between  aircraft 
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and  missile  areas.  The  Change  Engagement  Status  subroutine 
provides  the  necessary  program  bookkeeping  for  logical  consistency 
when  the  coordinated  missile  firing  procedure  is  used. 

The  dynamic  process  of  simulating  one  air  battle  from  start 
to  finish  forms  the  executive  routine  for  the  computer  simulation. 
This  executive  routine  consists  of  nothing  more  than  two  program 
subroutines  referred  to  as  SNE  and  TNE.   SNE,  Store  Next  Event,  is 
the  subroutine  that  takes  the  generated  information  pertaining  to 
an  interaction  and  properly  places  this  information  on  the  Event 
Store  List.  TNS,  Take  Next  Event,  is  the  subroutine  that,  at  the 
completion  of  any  of  the  five  events,  interrogates  the  information 
on  the  Event  Store  List  and  transfers  control  of  the  computer 
program  to  the  proper  subroutine. 

General  flow  charts  describing  the  logic  included  in  each 
event  of  the  simulation  plus  the  interrelationship  of  events  is 
included  as  Figure  4  through  Figure  9. 
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7 .   INPUT 

All  input  to  the  computer  program  is  contained  on  four 
punched  cards.  The  format  for  this  information  and  the  units 
to  be  used  are  detailed  in  Figures  10,  11,  and  12.  With  one  set 
of  input  cards  a  model  user  can  simulate  the  interactions  involved 
with  a  fixed  set  of  missile  areas,  with  fixed  values  describing 
the  missile  systems,  and  sets  of  aircraft  of  increasing  numbers. 
The  simulation  is  repeated  using  different  sets  of  aircraft  and 
different  sequences  of  random  numbers  in  order  to  provide  data  for 
statistical  analysis.   In  order  for  the  model  user  to  vary  parameter 
values  other  than  the  number  of  aircraft  a  new  set  of  input  cards 
must  be  supplied  to  the  computer  program.  With  the  proper  input 
values  the  model  can  be  made  deterministic  and  used  for  straight 
forward  calculations  of  sequences  of  intercept  points,  etc.  A 
complete  description  of  each  input .variable  follows  Figure  12. 

For  input  variables  whose  names  begin  with  I,  J,  K,  L,  M,  or 
N  the  input  values  must  be  integers  in  the  units  specified.  For 
all  other  input  variables  the  input  values  are  floating  point  and 
decimal  points  can  be  either  omitted  or  used  ss   desired  in  the 
space  allotted.  For  example  the  input  variable  name  AMAX  representing 
the  maximum  altitude  of  aircraft  is  a  floating  point  value  whose 
units  are  thousands  of  feet.   In  the  three  columns  allotted  for' 
AMAX  on  Card  3  the  model  user  could  enter  Blank-5-Q  or  5-0- .  for 
fifty  thousand  feet  and  could  also  enter  Blank-. -5  for  five  hundred 
feet. 
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COLUMNS 

NAME 

RANGE 

UNITS 

1 

NMISAR 

*  3 

■ 

2-3 

MTRF  ( 1) 

<  100 

4-5 

MTRF  (2) 

<  100 

6-7 

MTRF  (3) 

<  100 

8 

MAL(l) 

<  10 

9 

MAL(2) 

<  10 

10 

MAL(3) 

<  10 

11-13 

AX(1) 

<  1000 

MILES 

14-16 

AX(2) 

<  1000 

MILES 

17-19 

AX(3) 

<  1000 

MILES 

20-22 

AY(1> 

<  1000 

MILES 

23-25 

AY(2) 

<  1000 

MILES 

26-28 

AY(3) 

•<  1000 

MILES 

29-31 

RMR(l) 

<  1000 

MILES 

32-34 

RMR(2) 

<  1000 

MILES 

35-37 

RMR(3) 

<  1000 

MILES 

38-40 

RMAX(l) 

<  1000 

MILES 

41-43 

RMAX(2) 

<  1000 

MILES 

44-46 

RMAX(3) 

<  1000 

MILES 

47-50 

AVS(l) 

<  10000 

MILES /HR 

51-54 

AVS(2) 

<  10000 

MILES /HR 

55-58 

AVS(3) 

<  10000 

MILES /HR 

59 

IEND 

1  OR  0 

60-72 

NAME  OF 

PROGRAM  USER 

ALL  INPUT  INFORMATION  IS  RIGHT -JUSTIFIED  IN  THE  SPACE  ALLOTTED, 
FORMAT  FOR  INPUT  DECK  CARD  1 
FIGURE  10 
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3&- 


COLUMNS 

name 

RANGE 

UNITS 

1-2 

PK(D 

<  100 

% 

3-4 

PK(2) 

<  100 

% 

5-6 

PK(3) 

<  100 

z 

7-9 

ATM(1> 

<  1000 

MINUTES 

10-12 

ATM(2> 

<  1000 

MINUTES 

13-15 

ATM(3) 

<  1000 

MINUTES 

16-18 

ATX(l) 

<  1000 

MINUTES 

19-21 

ATX(2) 

<  1000 

MINUTES 

22-24 

ATX  (3) 

<  1000 

MINUTES 

25-27 

ASM(l) 

<  1000 

MINUTES 

28-30 

ASM(2) 

<  1000 

MINUTES 

31-33 

ASM(3) 

<  1000 

MINUTES 

34-36 

ASX(l) 

<  1000 

MINUTES 

37-39 

ASX(2) 

<  1000 

MINUTES 

40-42 

ASX(3) 

<  1000 

MINUTES 

43-45 

RTM(l) 

<  1000 

MINUTES 

46-48 

RTM(2) 

<  1000 

MINUTES 

49-51 

RTM(3) 

<  1000 

MINUTES 

52-54 

RTX(l) 

<  1000 

MINUTES 

55-57 

RTX(2) 

<  1000 

MINUTES 

58-60 

RTX(3) 

<  1000 

MINUTES 

ALL  INPUT  INFORMATION  IS  RIGHT -JUSTIFIED  IN  THE  SPACE  ALLOTTED. 

FORMAT  FOR  INPUT  CARD  2 
FIGURE  11 
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COLUMNS 

NAME 

RANGE 

UNITS 

1-3 

NBMAX 

£ 

20 

4-7 

BSX 

< 

10000 

MILES /HR 

8-11 

BSM 

< 

10000 

MILES /HR 

12-14 

BTX 

< 

1000 

MINUTES 

15-17 

BTM 

< 

1000 

MINUTES 

18-20 

NBOMB 

£ 

20 

21-23 

NBDEL 

<. 

20 

24-26 

GX 

< 

1000 

MILES 

27-29 

GY 

< 

1000 

MILES 

30-32 

RCC 

< 

1000 

MILES 

33-35 

CPA 

< 

180 

DEGREES 

36-37 

TMAX 

< 

100 

HOURS 

38 

OUTPUT 

1 

OR  0 

39-41 

NRl 

< 

1000 

42-44 

NRDEL1 

< 

1000 

45-47 

NR2 

< 

1000 

48-50 

NRDEL2 

< 

1000 

51-52 

NRPL 

£ 

20 

53-55 

AMAX 

< 

100 

K  FEET 

56-58 

AMIM 

< 

100 

K  FEET 

CARD  4 


COLUMNS 

NAME 

RANGE 

UNITS 

1-3 

STEP 

MINUTES 

4-5 

NFLAG 

0,1  OR  2 

6-7 

IHIST 

0   OR    1 

ALL  INPUT  INFORMATION  IS  RIGHT -ADJUSTED  IN  THE  SPACE  ALLOTTED, 
FORMAT  FOR  INPUT  DECK  CARD  3  AND  4 

FIGURE  12 

24 


DESCRIPTION  OF  INPUTS 


a 


CARD  1: 

NMISAR 

MTRF(I) 

MAL(I) 


AX(I) 

AY(I) 

RMR(I) 

RMAX(I) 

AVS(I) 

IEND 


NAME 


Number  of  missile  areas  in  the  game  simulation. 

Number  of  tracking  radars  assigned  to  missile 
area  I. 

Number  of  missile  launchers  assigned  to  missile 
area  I.  The  model  assumes  one  missile  launcher 
represents  the  salvo  size.  The  salvo  size  is 
not  explicitly  included  as  an  input  value. 

X-coordinate  of  missile  area  (I) . 

Y-coordinate  of  missile  area  (I) . 

Search  Radar  Maximum  Range  of  missile  area  (I) . 

Missile  Maximum  Range  for  missile  area  (I) . 

Average  Speed  of  missile  for  missile  area  (I) . 

To  change  the  values  of  any  input  parameters, 
except  the  number  of  aircraft,  another  set  of 
4  input  cards  is  required.  The  new  set  follows 
the  original  set  of  4  cards  in  the  data  deck 
ad-inf initum.  IEND  is  a  program  flag  to  indicate 
the  last  set  of  data. 

0  =>  Last  set  of  data 

1  =*  Another  set  of  4  input  cards  to  follow 

Programmer  name  or  identification.  Maximum  of 
13  Hollerith  characters. 


Note:   If  the  number  of  missile  areas  used  is  less  than  three  the 

computer  program  assumes  the  missile  areas  are  identified  in 
sequence  beginning  with  missile  area  1.  For  example,  if  only 
two  missile  areas  are  being  used  the  input  data  must  correspond 
to  missile  areas  1  and  2. 
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CARD  2: 
PK(I) 

ATM(I) 
ATX  (I) 
ASM  (I) 
ASX(I) 
RTM(I) 
RTX(I) 


Missile  salvo  Kill  Probability  for  missile  area  (I) . 
Salvo  size  is  not  explicitly  included  in  the  model. 
Salvo  size  is  included  by  coordinating  the  number 
of  missile  launchers,  MAL(I) ,  with  PK(I) . 

Minimum  Acquisition  Time  for  missile  area  (I)  .  -\ 

Maximum  Acquisition  Time  for  missile  area  (I). 

Minimum  Assessment  Time  for  missile  area  (I) . 

Maximum  Assessment  Time  for  missile  area  (I)  • 

Minimum  Reload  Time  for  missile  area  (I)  . 

Maximum  Reload  Time  for  missile  area  (I)  . 
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CARD  3: 

NBMAX 
BSX 

BSM 

BTX 
BTM 


I 


Maximum  Number  of  Aircraft  in  the  last  run. 

Maximum  Aircraft  Speed. 

Minimum  Aircraft  Speed.  For  all  aircraft  to 
enter  at  the  same  speed  BSX  and  BSM  should 
be  the  same  value. 

Maximum  Time  Spacing  between  aircraft. 

Minimum  Time  Spacing  between  aircraft.  For 
all  aircraft  to  enter  at  the  same  time  BTX  and 
BTM  should  be  set  to  zero.  For  all  aircraft 
to  enter  at  a  definitely  fixed  time  spacing 
between  aircraft  BTX  and  BTM  should  be  the 
same  value. 

NBOMB  Number  of  Aircraft  in  the  first  run. 

NBDEL  The  increment  used  to  increase  the  number  of 

aircraft  for  a  new  run.  The  number  of  runs  to 
be  made  with  one  set  of  input  cards  will  be  the 
minimum  integer  value  of  N  such  that 

NBOMB  +  (N)  (NBDEL)  ^  NBMAX 

GX  The  X-coordinate  of  the  center. 

GY  The  Y-coordinate  of  the  center. 

RCC  The  length  of  the  Radius  of  the  playing  area. 

CPA  The  Central  Playing  Angle  defining  the  playing 

area  and  aircraft  line  of  entry. 

TMAX  Maximum  Game  Time.  The  simulation  will  automatically 

end  when  no  more  interactions  are  possible.   If 
however  the  model  user  would  only  prefer  results 
for  a  fixed  period  of  time,  that  time  can  be 
entered  here.   If  complete  results  are  desired, 
use  a  very  large  number  such  as  9S.  See  note. 

OUTPUT  Program  flag  used  to  indicate  if  Standard  and/ 

or  Summary  output  is  desired. 

0:  Standard  and  Summary  output 
1:  Summary  only 

Note:  The  first  aircraft  enters  the  playing  area  at  time  =  1  hour. 
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CARD  3:  (Continued) 


OUTPUT 

(Continued) 


NRl 


NRDEL1 


NR2 


NRDEL2 


NRPL 
AMAX 
AMIN 


Normally  only  Summary  output  is  desired.  The 
use  of  Standard  output  greatly  increases  the 
computer  running  time. 

The  position  of  RNl,l,  the  first  random  number 
selected  for  the  simulation.  Since  this 
position  is  reached  by  generating  and  discarding 
the  numbers  in  previous  positions  use  a  small 
value,  say  200,  to  save  computer  time. 

The  increment  for  the  sequential  jump  on  NRl  for 
a  new  replication.  If  NRl  and  NRDELl  are  the 
same  on  different  sets  of  input  cards,  the  same 
sequences  of  random  numbers  will  be  used  for 
the  simulations  represented  by  the  different 
sets  of  input  cards. 

The  position  of  the  first  Random  Number  used  for 
the  generation  of  aircraft  entry  points,  times, 
speeds  and  altitudes. 

The  increment  for  the  sequential  jump  on  NR2  for 
a  new  replication.  If  NR2  and  NRDEL2  are  the 
same  on  different  sets  of  input  cards,  the  same 
sequences  of  random  numbers  will  be  used  to 
generate  the  sets  of  aircraft  for  the  simulations 
represented  by  the  different  sets  of  input  cards. 
If  it  is  the  'Case  that  all  other  aircraft  and' 
playing  area  information  is  the  same,  then  the 
sets  of  aircraft  used  in  the  simulations  will  be 
exactly  the  same. 

The  number  of  Replications  to  be  played. 

Maximum  Aircraft  Altitude. 


\J 


/ 


• 


Minimum  Aircraft  Altitude.  Aircraft  altitude  is 
used  only  to  compute  the  radar  horizon  using  the 
expression     ,— 

RH  =  1.25  "V  A    where 


RII  is  Radar  Horizon  in  miles  and 

A  is  Aircraft  altitude  in  feet 
An  altitude  of  zero  yields  a  radar  horizon  of  zero 
and  consequently  for  aircraft  with  an  altitude  of 
zero  no  missile  area  interactions  will  occur.  For 
all  aircraft  to  have  the  same  altitude  AMAX  and 
AMIN  should  be  the  same  value. 
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CARD  4: 


STEP 


NF1AG 


IHIST 


Step  is  a  time  step  value  used  for  computer 
efficiency.  When  a  missile  salvo  can  not  be 
fired  at  an  aircraft  because  of  (a)  no  launcher 
is  loaded,  (b)  no  tracking  radar  is  free,  or 
(c)  if  the  missile  area  fired  now  the  computed 
intercept  point  would  be  beyond  missile  maximum 

.  range,  then  a  new  Fire  Event  is  stored  at  current 
time  plus  the  value  of  Step.  This  value  must  not 
be  zero.  For  normal  aircraft  speeds  a  decent 
value  of  Step  is  probably  one  minute. 

This  value  controls  the  missile  firing  procedures: 

0:  Both  uncoordinated  and  coordinated 

procedures  are  used  in  the  simulation 

1:  Uncoordinated  procedure  used  only 
2:  Coordinated  procedure  used  only 

See  the  previous  section  on  DOCTRINE  for  a 

detailed  discussion. If  in  the  play  of  the  simulation 

only  one  missile  area  or  one  aircraft  is  used,  or  if  the 

spheres  of  influence  of  the  missile  areas  do  not  overlap 

then  the  uncoordinated  procedure  results  will 

be  the  same  as  the  coordinated  and  the  value 

of  NFIAG  should  be  1  in  order  to  save  computer 

time. 

This  value  controls  the  battle  history  output. 

0:  No  battle  history  output 

1:  A  battle  history  output  is  printed 

The  use  of  the  battle  history  output  greatly 
increases  computer  running  time  and  should  be. 
used  judiciously  if  at  all.  Therefore  this  value 
should  normally  be  set  to  zero. 
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8.   OUTPUT 

As  mentioned  earlier,  the  computer  program  presents  three 
different  sets  of  simulation  results;  Battle  History,  Standard 
and  Summary  output.  The  Battle  History  and  Standard  output  appear 
for  each  replication  and  are  printed  only  on  request,  i.e.,  when 
the  input  values  for  IHIST  and  OUTPUT  are  properly  set.   See  Section 
7  for  the  details  of  input.  The  Summary  output  is  always  printed 
for  each  run.  Computer  running  time  is  greatly  increased  with  the 
use  of  the  Battle  History  or  Standard  output.  This  section  presents 
the  format  of  the  output  so  that  the  model  user  may  interpret  the 
simulation  results. 

The  Battle  History  is  a  list  of  all  the  events  stored  by  SNE 
and  executed  by  TNE  during  the  replication  in  the  order  in  which 
the  events  were  stored  and  executed.  The  Battle  History  format  is 
illustrated  in  Figure  13.  Each  line  of  four  numbers  corresponds 
to  an  event  and  must  be  decoded.  The  events  being  executed  by  TNE 
are  represented  by  the  lines  beginning  on  the  left  while  the 
indented  lines  represent  events  that  are  being  stored  as  a  result 
of  the  event  being  executed.  The  decoding  format  is  as  follows: 


Time  of  Event 


Type  of  Event 


Aircraft  # 


Missile  Area  # 


The  time  appears  in  hours  and  the  aircraft  and  missile  areas  are 
identified  by  number.  The  type  of  event  is  decoded  as  follows: 
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NUMBER  EVENT 

1  Reload 

2  Set  Tracking  Radar  Free 

3  Intercept 

4  Change  Engagement  Status 

5  Fire 

As  an  example  the  Battle  History  presented  in  Figure  13  could 
be  read  as  follows.  Fire  events  were  stored  for  Aircraft  1  and 
missile  areas  1  and  3  at  1.205  and  1.171  hours  respectively  and 
for  aircraft  2  and  missile  areas  1  and  2  at  1.256  and  1.210  hours 
respectively.  It  must  be  pointed  out  that  this  is  the  order  in 
which  the  events  were  stored  and  is  not  the  order  in  which  the  events 
appear  on  the  Event  Store  List.  The  order  of  the  events  on  the 
Event  Store  List  is  chronological.  The  generation  of  these  four 
events  corresponds  to  the  program  logic  illustrated  in  Figure   4 . 
The  next  line  is  not  indented  and  represents  a  TNE,  i.e.,  at  1.171 
hours  missile  area  3  attempted  to  fire  at  aircraft  1.  As  a  result 
of  this  Fire  event  an  Intercept  event  is  stored  at  1.203  hours 
and  the  launcher  is  scheduled  to  be  reloaded  at  1.175  hours.  The 
execution  of  the  Reload  and  Intercept  events  appears  next  since 
they  were  the  two  earliest  events  on  the  Event  Store  List.  As  a 
result  of  the  Intercept  event  a  Fire  event  is  stored  at  1.255 
hours  indicating  the  salvo  did  not  kill  the  aircraft  and  Set  Tracking 
Radar  Free  and  Change  Engagement  Status  events  are  also  stored  at 
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1.255  hours.  The  history  of  the  battle  can  thus  be  read  from 
the  Battle  History.  Notice  that  at  1.255  hours  missile  area  3  is 
scheduled  to  fire  at  aircraft  1  and  no  events  are  stored  as  a 
result  of  the  scheduled  firing;  this  indicates  that  at  1.255  hours 
aircraft  1  was  either  dead,  beyond  the  point  of  closest  approach 
or  the  game  geometry  is  such  that  a  solution  to  the  intercept 
equation  does  not  exist.  It  can  be  seen  that  aircraft  1  was  killed 
at  1.233  hours  by  missile  area  1;  this  is  indicated  since  no  Fire 
event  is  stored  as  a  result  of  the  Intercept  event  at  1.233  hours. 
The  Standard  output  format  is  illustrated  in  Figure  W  along 
with  an  explanation  of  all  the  information  presented.  The  example 
of  Figure  14  can  be  seen  to  agree  with  the  Battle  History  of 
Figure  13.  The  first  row  of  this  Standard  output  contains  infor- 
mation pertaining  to  aircraft  1  and  can  be  read  as  follows;  in 
replication  5  using  the  uncoordinated  missile  firing  procedure 
aircraft  1  entered  the  simulation  at  1.000  hours  at  X  =  169.47 
miles  and  Y  »  571.44  miles  with  a  speed  of  1241.62  miles/hour  and 
at  an  altitude  of  5.68  thousand  feet;  missile  area  2  could  not  fire 
at  aircraft  1  and  the  earliest  time  of  fire  from  missile  areas. 1 
and  3  was  1.205  and  1.171  hours  respectively;  a  total  of  two 
salvos  was  fired  at  aircraft  1,  one  each  from  missile  area  1  and 
3  and  the  aircraft  was  killed  by  a  salvo  fired  by  missile  area  1. 
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Qr      X 


^Ta 


J 


TV 


— =9- 


1.171 


1.205 

V^l.171 

1.256 

1.210 

5  1 


1.175 
1.203 


3 


1.205 


1.209 
1.210 


1.214 
1.233 


1.234 


1.255 
1.255 
1.255 
1.256 
1.285 
285 
287 
287 
289 
289 


1.203 

1.175 
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1 
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1 
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1 
1 
1 
1 
2 
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3 
1 
3 
3 
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2 
1 
2 
4 
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4 
2 
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3 
3 
3 
1 
3 
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1 
1 
2 
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UNCOORDINATED  ATTACK  RESULTS    REP   5 


BMR  SMA1  SMA2  SMA3  TOTAL  A/K  XCOORD  YCOORD   BRH 


110  1  2         1    169.47    571.44    216.51  \  -* 

2  0  10  12    225.47      72.93    216.51 


ALTITUDE    SPEED      TIME       TlME(l)    TIME(2)    TIME<3) 


5.68    1241.62    1.000      1.205  .000         1.171 

5.68    1121.85    1.043      1.256         1.210  .000 


where: 

BMR 
SMA.1 


SMA2 
SMA3 
TOTAL 

A/K 


XCOORD 

YCOORD 

BRH 

ALTITUDE 

SPEED 

TIME 

TIME(l) 


TIMS  (2) 
TIME  (3) 


The  aircraft  number 

The  total  number  of  salvos   fired  at  this  aircraft 

by  missile  area   number   1 
Same  as  SMA.1  for  missile  area  number  2 
Same  as  SMA.1  for  missile  area  number  3 
The  total  number  of  salvos  fired  at  this  aircraft 

by  all  missile  areas 
The  number  of  the  missile  area  killing  the  aircraft. 

A/K  =  0   indicates   the  aircraft  was  not  killed  and 

penetrated   the  defenses. 
The  X-coordinate  of  aircraft  entry  point,    in  miles, 

into   the  playing  area. 
The  Y-coordinate  of  aircraft  entry  point,    in  miles, 

into   the  playing  area. 
Aircraft  radar  horizon  in  miles 
The  aircraft  altitude  in  thousands  of   feet 
The  aircraft   speed   in  miles/hour 
The  aircraft  entry  time   in  hours 
The  earliest  possible  time,    in  hours,    that  missile  area 

1  can  fire  at   the  aircraft  due   to  radar  horizon,    search 

radar  maximum  range  and  missile  maximum  range.     A  zero 

value   indicates   this  aircraft  can  not  be   fired  at  by 

this  missile  area. 
Same  as  TH1S(1)    for  missile  area  number  2 
Same  as  TIME(l)    for  missile  area  number  3 


STANDARD  OUTPUT  FORMAT 
FIGURE  /4 
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The  Summary  output  presented  at  the  end  of  each  run  is 
illustrated  in  Figure  15  along  with  an  explanation  of  the  format 
used.  This  summary  consists  of  the  total  number  of  aircraft 
kills  and  the  total  number  of  salvos  fired  by  missile  areas  for 
each  replication.  The  sample  mean,  variance  and  standard  deviation 
of  the  data  presented  is  also  included.  It  can  be  seen  that  the 
last  line  of  blocks  1  and  3  corresponds  to  a  summary  of  the 
information  contained  in  the  examples  of  the  Battle  History  and 
Standard  output. 
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KILLS 
NBOMB    REP    MAI:  MA2    MA3    TOTAL      MEAN 


1. 
1. 

. 
1. 


1. 
1. 


1. 
2. 
2. 
1. 
2. 


.33 
.67 
.67 
.33 
.67 


VAR 

.22 
.22 
.22 
.22 
.22 


SD 


Ulock    / 


MAI 


MEAN 

VAR 

SD 


.60 
.24 
.49 


MA2 

.60 
.24 
.49 


MA3 

.40 
.24 
.49 


TOTAL 

1.60 
.24 
.49 


3/ocfc    2 


MAI 
2. 

MA2 
2. 

MA3 
2. 

TOTAL 
6. 

MEAN 
2.00 

1. 
1. 
1. 

1. 
1. 

2. 

1. 
1. 

3. 
3. 
3. 

1.00 
1.00 
1.00 

1. 

1. 

1. 

3. 

1.00 

• 

MAI 

MA2 

MA3 

1 

..20 
.16 
.40 

1.00 
.40 
.63 

1.40 
.24 
.49 

SHOTS  FIRED 
VAR 


TOTAL 

3.60 
1.44 
1.20 


SD 


.00 

.00 

.67 

.82 

.00 

.00 

.00 

.00 

.00 

.00 

3/oc/c  3 


3kcK    4 


utmruL-ry  Duipsb  rorrnQX 
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where  for  block  1: 


NBOMB     The  number  of  aircraft  entered  in  this  replication 

REP       The  number  of  the  replication 

MAl       The  total  number  of  aircraft  killed  by  missile 

area  1  in  this  replication 
MA2       Same  as  MAl  for  missile  area  2 
MA3       Same  as  MAl  for  missile  area  3 
TOTAL      The  total  number  of  aircraft  killed  by  all  missile 

areas  in  this  replication 
MEAN    .   The  sample  mean  of  (MA.1, MA.2, MA 3) 
VAR       The  sample  variance  of  (MAl,MA2,MA3) 
SD        The  sample  standard  deviation  of  (MAl,MA2,MA3) 


for  block  2: 


MEAN 


VAR 
SD 


The  sample  mean  of  the  number  of  aircraft  killed 
by  each  and  by  all  missile  areas.  The  entry  for 
MEAN  under  MAl,  5.80,  is  the  sample  mean  of  the 
column  of  numbers  under  MAl  in  block  1 
The  sample  variance  for  the  same  set  of  data 
The  sample  standard  deviation  for  the  same  set  of 
data 


Probability  of  Survival  ■ 


Number  of  replications  in  which 
all  aircraft  are  killed 

Total  number  of  replications 


and  represents  the  probability  that  an  aircraft 
attack  will  not  penetrate  the  missile  defenses. 

The  same  information  is  then  presented  for  the  total  number  of 
salvos  fired,  i.e.,  blocks  3  and  4. 


SUMMARY  OUTPUT  FORMAT 
FIGURE  15  (Continued) 
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9 .   REMARKS 

No  significant  results  have  been  obtained  to  date  with  the 
model  presented.  This  report  is  the  culmination  of  effort  to 
build  the  model  and  make  it  available  for  use  to  interested 
parties;  in  particular  to  the  students  of  Operations  Analysis 
at  the  U.  S.  Naval  Postgraduate  School.  A  parametric  study 
using  the  model  is  now  underway. 

The  appendices  to  this  report  contain  the  detailed  flow 
charts  of  the  computer  program,  the  CDC -FORTRAN -60  program  listing 
and  the  format  for  using  the  program  at  the  Postgraduate  School. 

The  author  thanks  Miss  Patricia  Hoang  for  her  able  assistance 
in  programming  and  debugging  the  computer  program  and 
Mr.  Cloyd  Marvin  of  The  Johns  Hopkins  University  Applied  Physics 
Laboratory  for  the  IBM- 7090  FAP  version  of  the  random  number 
generator. 
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APPENDIX  I 
Detailed  Flow  Charts 

This  appendix  contains  the  detailed  flow  charts  for  the 
computer  program.  The  flow  chart  pages  are  sequenced  from  A  to 
AG  and  the  flow  chart  symbols  for  logical  continuity  from  page 
to  page  are: 


and 

(70^S)-«1 

The  flow  chart   symbol  for   logical  continuity  on  the  same  page   is: 

►<!£> 
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APPENDIX  II 
Program  Names 

This  appendix  contains  the  name  and  definition  of  the  symbols 
appearing  in  the  computer  program  and  the  flow  charts.  All 
symbols  of  importance  are  included  except  those  already  defined 
in  Section  7.  The  names  in  this  appendix  are  listed  in  alphabetical 
order. 
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NAME 
ABAR(I,K) 

ALT(I,N) 
AM(I,J,K) 

AMT(I,K) 

APK(J) 
ASD(I,K) 

AVAR(I,K) 

BDPC(I,J,N) 

BM(J,K) 

BME(I,J) 


BMT(K) 

BRH(I,N) 

BS(I,N> 

BSD(J,K) 

BSDT(K) 


DEFINITION 

The  sample  mean  for  the  number  of  kills  by  all 
missile  areas  in  replication  I.* 

The  altitude  of  aircraft  I  in  replication  N. 

The  number  of  kills  for  missile  area  J  in 
replication  I.* 

The  total  number  of  kills  for  all  missile  areas 
in  replication  I.* 

The  salvo  kill  probability  for  missile  area  J. 

The  sample  standard  deviation  for  the  number  of 
kills  by  all  missile  areas  in  replication  I.* 

The  sample  variance  for  the  number  of  kills  by  all 
missile  areas  in  replication  I.* 

The  distance  from  missile  area  J  to  the  point  of 
closest  approach  for  aircraft  I  in  replication  N. 

The  sample  mean  for  the  number  of  kills  by  missile 
area  J  in  one  run.* 

A  non-zero  value  indicates  that  aircraft  I  is 
currently  engaged  by  missile  area  J.  Before 
intercept  the  value  is  the  time  of  intercept  and 
after  intercept  the  value  is  the  time  of  intercept 
plus  assessment  time.  A  zero  value  indicates  no 
current  engagement. 

The  sample  mean  for  the  number  of  kills  for  all 
missile  areas  for  one  run. 

The  radar  horizon  for  aircraft  I  in  replication  N. 

The  speed  of  aircraft  I  in  replication  N. 

The  sample  standard  deviation  for  the  total  number 
of  kills  for  missile  area  J  in  one  run.* 

The  sample  standard  deviation  for  the  total  number 
of  kills  for  all  missile  areas  in  one  run.* 
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/* 


NAME 
BT(I,N) 

BTMA(I,J,N) 

BT?C(I,J,N) 

BV(J,K) 

BVT(K) 

BVX(I,N) 

BVY(I,N) 

BX(I,N) 

BY(I,N) 

CBAR(I,K) 

CSD(I,K) 

CTOT(I,K) 

CTOTS(I,J,K) 

CVAR(I,K) 

DM(J,K) 

DMT(K) 


DEFINITION 

The   tirna  aircraft  I  enters   the  playing  area   in 

replication  N. 

The  earliest  possible  time  missile  area  J  can  fire 
at  aircraft  I  in  replication  N. 

The  time  aircraft  I  is  at  the  point  of  closest  approach 
for  missile  area  J  in  replication  N. 

The  sample  variance  for  the  total  number  of  kills  for 
missile  area  J  in  one  run.* 

The  sample  variance  for  the  total  number  of  kills 
for  all  missile  areas  in  one  run.* 

The  X-coraponent  of  velocity  for  aircraft  I  in 
replication  N. 

The  Y-cornponent  of  velocity  for  aircraft  I  in 
replication  N. 

The  X-coordinata  of  the  entry  point  into  the  playing 
area  for  aircraft  I  in  replication  N. 

The  Y-coordinate  of  the  entry  point  into  the  playing 
area  for  aircraft  I  in  replication  N. 

The  sample  variance  for  the  number  of  salvos  fired "by  all 
missile  areas  in  replication  I.* 

The  sample  standard  deviation  for  the  number  of  salvos 
fired  by  all  missile  areas  in  replication  I.* 

The  total  number  of  salvos  fired  by  all  missile  areas 
in  replication  I.* 

The  total  number  of  salvos  fired  by  missile  area  J  in 
replication  I.* 

The  sample  variance  for  the  number  of  salvos  fired  by 
all  missile  areas  in  replication  I.* 

The  sample  mean  for  the  number  of  salvos  fired  by  missile 
area  J  in  one  run.* 

The  sample  mean  for  the  number  of  salvos  fired  by  all 
missile  areas  in  one  run.* 
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NAME 
DSD(J,K) 

DSDT(K) 

DV(J,K) 

DVT(K) 

EVS(I) 
IBA(I) 

IBMEC(I,J) 

IEVS1(I) 

IEVS2(I) 

IEVS3(I) 

IE1 
IE2 
IE3 

IFE(I) 

IKILL(I) 

IRUN 


DESCRIPTION 

The  sample  standard  deviation  for  the  number  of  salvos 
fired  by  missile  area  J  in  one  run.* 

The  sample  standard  deviation  for  the  total  number  of 
salvos  fired  by  all  missile  areas  in  one  run.* 

The  sample  variance  for  the  number  of  salvos  fired  by 
missile  area  J  in  one  run.* 

The  sample  variance  for  the  total  number  of  salvos 
fired  by  all  missile  areas  in  one  run.* 

The  time  of  the  Ith.  event  in  the  Event  Store  List. 

A  non-zero  value  indicates  that  aircraft  I  has  been 
killed  in  the  current  replication. 

The  number  of  salvos  fired  at  aircraft  I  by  missile 
area  J  in  the  current  replication. 

The  type  of  event  in  the  Ith.  position  in  the  Event 
Store  List. 

The  aircraft  number  associated  with  the  Ith.  event  in 
the  Event  Store  List. 

The  missile  area  number  associated  with  t he  Ith.  event 
in  the  Event  Store  List. 

The  type  of  the  current  event. 

The  aircraft  number  associated  with  the  current  event. 

The  missile  area  number  associated  with  the  current 
event. 

The  total  number  of  initial  FIRE  events  stored  in 
replication  I. 

A  non-zero  value  indicates  the  number  of  the  missile 
area  that  killed  aircraft  number  I  in  the  current 
replication. 

A  0  or  1  value  indicates  respectively  that  the  program 
is  executing  an  uncoordinated  or  coordinated  replication. 
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NAME 

ISl 

IS2 

IS3 

ITOT(I) 

MAL(J) 

MALO(J) 

MTRF(J) 

MTRFO(J) 

NOL 

NREP 

NRlO 
NR20 
NX 
PS(K) 

RN 

RNG(A) 


DESCRIPTION 

The  type  of  event  being  stored. 

The  aircraft  number  associated  with  the  event  being 
stored. 

The  missile  area  number  associated  with  the  event 
being  stored. 

The  total  number  of  salvos  fired  by  all  missile 
areas  at  aircraft  I  in  the  current  replication. 

The  current  number  of  loaded  missile  launchers  at 
missile  area  J 

The  original  number  of  loaded  missile  launchers  at 
missile  area  J. 

The  current  number  of  available  tracking  radars  at 
missile  area  J. 

The  original  number  of  available  tracking  radars  at 
missile  area  J. 

The  number  of  events  currently  on  the  Event  Store  List. 

A  non-zero  value  indicates  a  new  set  of  aircraft  will 
be  generated  for  the  next  execution  of  the  program. 

The  original  value  of  NRl. 

The  original  value  of  NR2. 

The  usual  argument  for  RNG(NX)  with  a  value  of  1. 

The  probability  of  survival  for  the  Summary  Output. 

The  functional  value  of  RNG(A) ,  i.e.,  the  value  of 
the  computer  generated  random  number. 

The  random  number  generating  function.  The  value  of  ' 
A  determines: 

A  =  0  :   generate  the  very  first  number  in  the 

sequence. 
A  =  k  ^  0:   generate  the  number  k  positions  from 
the  last  number  generated. 
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NAME  DESCRIPTION 

Si  The  time  of  the  event  being  stored. 

T  The  time  of  the  current  event. 

TF  The  tirae  of  flight  for  the  missile  as  calculated 

in  the  FIRE  event. 


NOTE:  *:  K  ■  0  or  1  indicates  that  the  replication  or  run  being 
processed  is  respectively  uncoordinated  or  coordinated. 
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APPENDIX  III 
Computer  Program 

This  appendix  contains  the  complete  computer  program.  The 
program  is  in  FORTRAN-60  for  the  CDC-1604.  All  numerical  restrictions 
and  variable  names  are  such  that  the  program  is  compatible  with  IBM 
FORTRAN-II  for  the  IBM-70S0-94.   On  the  CDC-1604  the  program  does 
not  compile  in  one  pass  because  of  space  allocations. 

The  random  number  generator,  FUNCTION  RNG(NR) ,  appearing  at  the 
end  of  the  program  is  a  combination  of  both  the  assembly  and  compiler 
language  for  the  CDC-1604.  This  function  would  have  to  be  reprogrammed 
for  the  IBM-7090-94.  The  logic  description  of  this  function  is  as 
follows: 

RNG  is  a  random  number  generator  subroutine  to  be  us.ed  with 
FORTRAN  on  the  IBM  7090  or  IBM  7094.   It  produces  a  floating  point 
random  number  which  is  uniformly  distributed  over  the  open  interval 

(0,  1).  The  generator  produces  random  numbers  according  to  the  cycle 

15      35 
X   -  X   x  5   (MOD  2  > 
i      i-1 

RN  -  2~35  x  X 
i 

15 
where  X   =  5 

o 

The  routine  is  a  standard  FORTRAN  function  -  type  subprogram,  . 

named  RNG (MR) .   Its  argument  I  is  a  fixed  point  variable  or  fixed 

point  constant  which  tells  how  many  times  to  cycle  the  generator 

from  its  previous  value.   If  I  =  0,  the  generator  is  reset 
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15v  15    -35 

(X  is  set  to  5  )  and  the  first  number  is  5   x  2 
o 

The  example  below  is  given  as  RNG(NR)  might  be  used. 

1  A  -  RNG(O) 

2  B  -  RNG(50) 

3  C  =  RNG(l) 

4  D  =  RNG(l) 

15   -35 

The  first  statement  causes  the  generator  to  be  reset  and  5  x  2 

is  put  in  A.   In  the  second  statement,  the  generator  is  cycled  fifty 
times  from  its  previous  value  and  the  floating  point  result  is 
placed  in  B.  The  last  two  statements  place  the  next  two  numbers 
generated  in  C  and  D  respectively. 
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PROGRAM  REGAME 

DIMENSION  IBA(20) , IBMEC ( 20,3 ) ,BME ( 20,3 ) ,BVX( 20,20 ) » 

1  BVY(20»20)  ,APK(3)  »MTRF(  3  )  »  I  KILL(  20  )  » ITOT(  20)  » 

2  MTRFO(3) ,MALO(3) »AM( 20,3,2) »AMT(20,2), 

3  ABAR(20,2) »AVAR(20>2) »ASD(20,2) »BM(3.2) , 

4  BMT(2)»BV(3»2)»BVT(2)»BSD(3»2)  »BSDT (2 )  ,PS ( 2  )  » 

5  EVS(200) ,IEVS1(200),IEVS2(200) »IEVS3(200> » 

6  AX(3) ,AY(3) »MAL(3) »RMR(3) ,RMAX(3) ,AVS(3)» 

7  ATM(3),ATX(3) ,ASM(3) ,ASX(3) *RTM(3) 
DIMENSION  RTX(3),PK(3) »BX ( 20 »20 ) »BY( 20,20 ) »BS ( 20»20 ) , 

1  BT(20,20),BRH(20,20) ♦ ALT ( 20*20 ) ,BTMA( 20 ,3,20) ♦ 

2  BTPC(20,3,20) ,BDPC ( 20 ,3,20 ) »CTOTS( 20, 3,2 ) , 

3  CTOT(20,2) ,CBAR(20,2) ,CVAR(20,2) ,CSD(20»2) • 

4  DM (3,2) ,DMT(2),DV(3,2) ,DVT( 2 ) ,DSD( 3,2 ) • 

5  DSDT(2) ,NAME(2) ,IFE(20) 
C 

C      READ  INPUT 
C 

1  READ  INPUT  TAPE  5,2    »NMI SAR , ( MTRF ( I ) , I =1 ,3 ) , (MAL( I ) » 

1  1=1,3) »(AX( I) , 1=1,3), (AY ( I) ,I=1»3)»(RMR(I) ,1*1,3) 

2  ,(RMAX( I ) ,1=1,3  )»(AVS( I ) ,1  =  1,3) , I  END, NAME ( 1 ) , NAME( 2 ) 

2  FORMAT  (11,312, 311,  12F3.0 , 3F4.0 ♦ 1 1 »A8, A5 ) 

READ  INPUT  TAPE  5,3    , (PK ( I ) , 1=1 ,3 ) • ( ATM( I ) , 1=1 ,3) , 

1  (ATX(I),I=1,3) ,(ASM(I),I=1,3),(ASX(I) ,1=1,3), 

2  (RTM( I ) ,1=1,3) ,(RTX(I ),I=1,3) 

3  FORMAT  (3F2.2,18F3.0) 

READ  INPUT  TAPE  5,4    ,NBMAX »BSX,BSM,BTX ,BTM,NBOMB, 

1  NBDEL, GX,GY,RCC,CPA,TMAX, OUTPUT, NR1,NRDEL1, 

2  NR2,NRDEL2,NRPL,AMAX,AMIN 

4  FORMAT  U3»2F4.0,2F3.0»2I3,4F3.0»F2.0»I1»4I3»I2,3P2F3»0) 
READ  INPUT  TAPE  5,  494,  STEP*  NFLAG,  IHIST 

494  FORMAT  (F3.0,  212) 

PRINT  9991,  NAME(1)»NAME(2) 


9991 

FORMAT  (3X.2A8) 

C 

C 

SET  CONSTANTS 

c 

IF  (NFLAG-1) 

495.495, 

495 

IRUN=0 
GO  TO  497 

496 

IRUN=1 

497 

i 

NRRN=0 

N=l 

NREP=0 

NOL   =  0 

DO  5  1=1,20 

IKILL(I)  =0 

IBA(I)    =0 

ITOT(I)   =0 

IFE(I)=0 

DO  5  J«lt3 

BME(I,J)=0. 

IBMEC(I»J)=0 

5 

CONTINUE 

IF(NMISAR-3) 

6,8,9999 

496 
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8 
9 

10 


L=NMISAR+1 

DO  7  K»lt2 

DO  7  J  =  Lt3  • 

BM( JtK)=0. 

BV(JtK)=0. 

BSD(JtK)=0. 

DM(JtK)«0. 

DV<JtK)=0. 

DSD(JtK)«0. 

DO  7  I«ltNRPL 

AM( ItJtK)=0, 

CTOTSUtJtK)«0. 

CONTINUE 

PRINT  INPUT 

WRITE  OUTPUT  TAPE  6t  9    tNAME ( 1 ) tNAME( 2 ) 
FORMAT  U1H1INPUT  DATA »75X ,2A8/// ) 
WRITE  OUTPUT  TAPE  6,10   t  GX tGYtRCCtCPA 


FORMAT ( 5X » 15HGAME  PARAMETERS// » 

1  30H  XCOORDINATE  OF  CENTER 

2  30H  YCOORDINATE  OF  CENTER 

3  31H  RADIUS  OF  CIRCLE 

4  31H  CENTRAL  PLAYING  ANGLE 


lOXt 

(GX) 

(GY) 

(RCC) 

(CPA) 


F8.2, 

F8.2t 

■  F7.2 

'  F7.2 


2X»5HMILES/10X. 
2X»5HMILES/10X» 
t2Xt5HMILES/10Xt 
t2Xt  7HDEGREES////) 


11 

12 

13 
14 
15 
16 

17 
18 

19 

20 

21 

22 


1 


WRITE 
FORMAT 

WRITE 

FORMAT 

WRITE 

FORMAT 

WRITE 

FORMAT 

WRITE 

FORMAT 

WRITE 

FORMAT 

ANGE) 

WRITE 

FORMAT 

WRITE 

FORMAT 

) 

WRITE 

FORMAT 

NIMUM) 

WRITE 

FORMAT 

XIMUM) 

WRITE 

FORMAT 

IMUM) 

WRITE 

FORMAT 

IMUM) 


OUTPUT 
(25H  MI 
8Xt  1 
OUTPUT 

(21X. 
OUTPUT 
(21Xt  3 
OUTPUT 
<21Xt  I 
OUTPUT 

<21XtI 
OUTPUT 
(21Xt  3 

OUTPUT 
(21Xt  3 
OUTPUT 
(20X*  3 

OUTPUT 
(21X,  3 

OUTPUT 
(21X,  3 

OUTPUT 
(21X*  3 

OUTPUT 
(21X.  3 


TAPE  6, 
SSILE  A 
Hit  9Xt 
TAPE  6t 
3(F6,2t 
TAPE  6t 
(F6»2t4 
TAPE  6t 
3t7Xt  I 
TAPE  6, 
3t7XtI3 
TAPE  6t 
(F6.2t4 


11 
REA  INFORMATION//lXt  16H 

lH2t  9Xt  1H3/) 

12   t  (AX( I) t  I»lt3) 
4X)t  5HMILESt  7Xtl9HAX 


13   t  (AY(I 
X)t  5HMILES 

14  t  (MAL( 
3t  7Xt  I3t 

15  t  (MTRF 
t7XtI3tl9Xt 

16  t  (RMR( 
X)t  5HMILES 


TAPE  6t  17  ♦  (RMAX 
(F6.2»4X) t  5HMILES 
TAPE  6,  18  t  (AVS( 
(F7.2t3X) tlXt  5HMI 

TAPE  6t  19  t  (ATM( 
(F6.2t4X)t  5HMINS 

TAPE  6t  20  t  (ATX( 
(F6.2t4X)t  5HMINS 

TAPE  6t21  t  (ASM( 
(F6.2t4X)t  5HMINS 

TAPE  6t  22  t  (ASX( 
(F6,2t4X)t  5HMINS 
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tl*lt3) 
7X»19HAY 

)  tl»lt3) 

9Xt   23HMAL 

I) tl=lt3) 
22HMTRF   =  NR 

)tl=lt3) 
7Xt  30HRMR 

I) »I=lt3) 

7X»25HRMAX 
) t I=lt3) 
/HRt7Xt25HAVS 

) tl=lt3) 
7Xt32HATM 

) tl»lt3) 
7Xt32HATX 

) tl=lt3) 
7Xt31HASM 

) tl=lt3) 
7Xt31HASX 


MISSILE  AREA  NR 

»  XCOORDINATE) 
»  YCOORDINATE) 

*  NR  OF  LAUNCHERS) 
OF  T  RADARS) 

*  SEARCH  RADAR  MAX  R 

=  MISSILE  MAX  RANGE) 
*  MISSILE  AVG  SPEED 

*  ACQUISITION  TIME  MI 
=  ACQUISITION  TIME  MA 
■  ASSESSMENT  TIME  MIN 

*  ASSESSMENT  TIME  MAX 


WRITE  OUTPUT  TAPE  6,   23  ♦  (RTM(I).  I=l»3) 
23  FORMAT(21X»  3(F6.2»4X)>  5HMINS  »  7X,27HRTM 
1) 
WRITE  OUTPUT  TAPE  6,   24  *  ( RTX < I ) ♦ I =1 t 3 ) 


=  RELOAD  TIME  MINIMUM 


3(F6.2»4X)»  5HMINS  ♦  7X»27HRTX    =  RELOAD  TIME  MAXIMUM 


t  (APK(I).  1=1.3) 
5HPER,  ♦  7X.34HPK 


24  FORMAT(21X» 
1) 

DO   25   I  =  l»3 

25  APK(D=100.*  PK(D 
WRITE  OUTPUT  TAPE  6>  26 

26  FORMAT(21X»  3<F6.2t4X). 
1R  SALVO) 

WRITE  OUTPUT  TAPE  6,  27 

27  FORMAT*///  18HBOMBER  INFORMATION//) 
WRITE  OUTPUT  TAPE  6» 

28  FORMAT(5X»21H  MAXIMUM 

1  5X.21H  MINIMUM 

2  5X»26H  MAXIMUM 

3  5X»26H  MINIMUM 

4  5X»31H  MAXIMUM 


=  KILL  PROBABILITY  PE 


5X.31H  MINIMUM 


28.  »  BSX»BSM,AMAX,AMIN,BTX»BTM 
SPEED  (BSX)=  F6.1»2X»  8HMILES/HR/ 
SPEED  (BSM)=  F6.1»2X»  8HMILES/HR/  ' 
ALTITUDE  (AMAX)=  F8. l»2Xt4HFEET/ 
ALTITUDE  (AMIN)=  F8. 1 t 2X.4HFEET/ 
SEPARATION  TIME  (BTX)=  F7.2»2X,  5HMINS./ 
SEPARATION  TIME  (BTM)=  F7. 2 »2X»5HMI NS./1H1 ) 


INITIALIZE 

DO  29  I*1»NMISAR 
ATM( I )=ATM( I )/60.0 
ATXU)=ATX( D/60.0 
ASM(D=ASM(  D/60.0 
ASXU  )=ASX(  D/60.0 
RTM(I)=RTM(  D/60.0 

29  RTX(I)»RTX(  D/60.0 
CPA=CPA/57. 29577 
BTX=BTX/60.0 
BTM=BTM/60.0 
STEP=STEP/60.0 
NX»1 

DO    30    J  =  l,3 
MALO(J)=MAL(J) 

30  MTRFOi J)»MTRF(J) 
NR10=NR1 
NR20=NR2 

COMPUTE  BOMBER  ENTRY  POINTS 


31  RN=RNG(0) 
RN=RNG(NR2) 

32  1  =  1 

33  RN«RNG(NX) 

IF  (RN-0.5)  34»36»35 

34  B=CPA*(0.5-RN) 
X=RCC*COSF(B) 
Y=RCC*SINF(B) 
BX(I»N)=GX+X 
BY(I»N)»GY+Y 
K=l 

GO  TO  37 

35  B»CPA*< RN-0.5) 
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X=RCC*COSF(B) 

Y=RCC*SINF(B) 

BX(  I»N)=GX+X 

BY( I,N)*GY-Y 

K  =  2 

GO  TO  37 

36  BX(I»N)=GX-fRCC 
BY(I»N)«GY 

K  =  3 

37  RN=RNG(NX) 
C 

C      COMPUTE  BOMBER  SPEED  AT  ENTRY 
C 

BS( I»N)=(BSX-BSM)*RN+BSM 

IF(K-1)39»38»39 

38  BVX(I.N)»-BS( I»N)*COSF(B) 
BVY< I»N)=-BS( I.N)*SINF(B) 
GO  TO  42 

39  IF  (K-2)  41,40»41 

40  BVXU»N)=-BS(  WN)*COSF(B) 
BVY(I»N)=BS(I»N)*SINF(B) 
GO  TO  42 

41  BVX(I,N)=-BS< I,N) 
BVY(ItN)=0. 

C 

C      COMPUTE  BOMBER  TIME  AT  ENTRY 

C 

42  IF(I-l)  44«43>44 

43  BT(I»N)=1.0 
GO  TO  45 

44  RN=RNG(NX)  ? 
BT(I»N)=(BTX-BTM)*RN+BTM+BT(I-1»N) 

C 

C      COMPUTE  BOMBER  ALTITUDE  AND  BOMBER  RADAR  HORIZON 

C 

45  RN=RNG(NX) 
ALT(I»N)=(AMAX-AMIN)«RN+AMIN 
BRH(I»N)=1.25*SQRTF(ALT(I»N) ) 
ALT(I»N)*ALT( I»N)/5280.0 

46  J  =  l 

47  A2=(BX( I»N)-AX< J) )**2+(BY< I  >N ) -AY< J ) ) **2 
A=SQRTF(A2) 

B2=(AX( J)-GX)**2  +  ( AY(J)-GY)**2 

B=SORTF(B2) 

C2=  RCC*RCC 

C  =  RCC 

CB=  <-B2+A2+C2)/<2.*A*C> 

SB2=  1«0-CB**2 

SB*SQRTF(SB2) 
C 

C      COMPUTE  DISTANCE  TO  PT  OF  CLOSEST  APPROACH  FROM  MISSILE  AREA 
C 

BDPC( I»J»N)=A*SB 

IFtBDPCU  #J»N)-RMAX( J)  > 48 • 48 ,57 

48  IF(A-RMRtJ))  55»55»49 

49  IFtBDPCU  •  J. N)-RMR  (J)  )  50*50,57 
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50  IF(BRH( I.N)-RMR(J) )505»51.51 

51  IF  (RMR(J)-RMAX(J) )  52,52,507 

52  Z2=RMR( J)**2-BDPC( I,J,N)**2 
GO  TO  53 

50  5  IF(BRH( I,N)-BDPC( I»J,N) ) 57, 506, 506 
7    506  IF  <BRH(I,N)-RMAX( J) )  508,508,507 

507  X=A*CB 
BTPC(I,J,N)=  X/BS( I,N)+BT( I,N) 
Z=BDPC(I»J,N)**2 
ZR2=RMR(J)**2  -  Z 
ZM2*  RMAX(J)**2  -  Z        • 
ZB2=  BRH( I,N)**2  -  Z 
ZM  *  SQRTF(ZM2) 
IF  (ZR2  -  ZB2)  800,  801,  801 

800  ZR=  SQRTF(ZR2) 
60  TO  802 

801  ZR=  SQRTF(ZB2) 

802  TRMAX=  (X-ZM) /BS( I ,N )  +  BT(I,N) 
TZR«  (X-ZR)/BS( I,N)  +  BT(I,N) 
Z*  RMAX(J)**2  +  ALT(I,N)**2 
Sl=  TRMAX  -  SQRTF(Z)/AVS(J) 

T  RN=RNG(NX) 

Z  =  (ATX( J)-ATM(J) )*RN  +  ATM(J) 
IF  ((Sl-Z)-TZR)  804,803,803 

803  BTMA(I,J,N)  ■  SI 
GO  TO  56 

804  Sl«  TZR+Z 
GO  TO  803 

508  Z2=BRH(I»N)**2-BDPC( I»J»N)**2 

53  Z=SQRTF(Z2) 
X=  A*CB 

C 

C      COMPUTE  TIME  AT  PT  OF  CLOSEST  APPROACH 

C 

BTPC(I,J,N)=X/BS( I»N)+BT(I ,N) 

Y  =  X-Z 
C 

C      TIME  MISSILE  AREA  ENTRY 
C 

BTMAU  ,J,N)=BT(  I ,N )+Y/BS( I ,N ) 

54  RN*  RNG(NX) 
S1=*BTMA( I,J,N)+(ATX( J)-ATM(J) )*RN+ATM(J) 
BTMA(I,J,N)=S1 
GO  TO  56 

55  IF  (A-BRH(I,N))  555,555,49 
555  BTMA  ( I»J,N)  =  BT ( I,N) 

X=SQRTF(A**2-BDPC( I,J,N)**2) 
BTPCUtJ»N)=X/BS(I»N)+BT( I,N) 
GO  TO  54 

56  IS1=5 
IS2«I 
IS3=J 

IFE(N)»IFE(N)+1 
ASSIGN  58  TO  NEXT 
GO  TO  700 

57  BTMAU  »J,N)»0. 
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58  IF(J-NMISAR)  59»60»9999 

59  J=J+1 
GO  TO  47 

60  IF(I-NBOMB)  61»62,9999 

61  1=1+1 
GO  TO  33 

62  RN=RNG(0) 
RN=RNG(NR1) 
GO  TO  600 

C 

C      TNE 

C 

600  IF(NOL)  9999.6000»601 

601  T=EVS(1) 
IE1=IEVS1(1) 
IE2*IEVS2<1) 
IE3=IEVS3(1) 

IF  (IHIST)  9999,  952»  950 

950  WRITE  OUTPUT  TAPE  6.  951.  T»  IE1»  IE2.  IE3 

951  FORMAT  (10X.  F7.3.  315) 

952  IF  (NOL-1)  9999*  604,  602 

602  DO  603  I=2tNOL 
EVS(I-l)  =  EVS(I) 
IEVSKI-1)  ■  IEVSKI) 
IEVS2U-1)  =  IEVS2U) 

603  IEVS3( 1-1)  =  IEVS3U ) 

604  EVS(NOL)  =0, 
IEVSKNOL)  =  0 
IEVS2(NOL)  =  0 
IEVS3(NOL)  =  0 
NOL  =  NOL  -  1 

IF  (  T  -  TMAX  )   605»605»608 

605  IF  (IE1  )  9999»9999t606 

606  IF  (IE1  -  5)  ,607t607,9999 

607  GO  TO  ( 1000»2000»3000»4000,5000)»IE1 

608  WRITE  OUTPUT  TAPE  6.609»  NBOMB.NRPLI 

609  FORMAT  (26H1MAXIMUM  TIME  EXCEEDED  FOR  I5t  8HBOMBERS  » 
1  4HRUN  13) 

GO  TO  6000 
C 

C      SNE 
C 

700  IF  (NOL)  9999»  701.  702 

701  M  =  1 

GO  TO  706 

702  M   =  NOL 

703  IF(S1  -  EVS(M))  704,709»708 

704  EVS(M+1)  =  EVS(M) 
IEVSKM+1)  =  IEVSKM) 
IEVS2(M+1)  =  IEVS2(M) 
IEVS3(M+1)  =  IEVS3(M) 
IF  (  M-l  )9999»706t705 

705  M  *  M-l 
GO  TO  703 

706  IF  (IHIST)  9999»  955»  953 

953  WRITE  OUTPUT  TAPE  6.  954»  Sit  ISlt  IS2t  IS3 
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954  FORMAT     (20Xf    F7.3»    315) 

955  EVS<M)=S1 
IEVSKM)    =    IS1 
IEVS2(M)    ■    IS2 
IEVS3(M)    *    IS3 
NOL    ■    NOL    +    1 

707  GO    TO   NEXT 

708  M    »    M    +    1 
GO    TO    706 

709  IFUS1    -    IEVSKM))       704»710»708 

710  IFUS2    -    IEVS2(M))       704»711»708 

711  IF(IS3    -    IEVS3(M))       704.707»708 
C 

C  RELOAD 

C 

1000    J    ■    IE3 

MAL(J)     ■    MAL(J)    +    1 

GO    TO    600 
C 

C  T    R    FREE 

C 
2000    J    =    IE3 

MTRF(J)    =    MTRF(J)    +    1 

GO    TO   600 
C 

C  INTERCEPT 

C 

3000  I=IE2 
J«IE3 
IF    (IBA(I))    600»3001»600 

3001  RN=RNG(NX) 
IF(RN-PK(J))    3002»3002t3009 

3002  IF(IRUN)    3013»3012»3013 
3012    K»l 

3003  IF(BME( I»K) ) 3004.3006 t 3004 

3004  RN=RNG(NX) 
S1=T+(ASX(K)-ASM(K))    *RN+ASM(K) 
BME(I»K)«S1 
IS3  =  K 
IS1»2 
IS2  =  I 

ASSIGN    3005    TO   NEXT 
GO    TO    700 

3005  IS1«=4 
IS2*I 

ASSIGN  3006  TO  NEXT 
GO  TO  700 

3006  IF(K-NMISAR)3007»3008»3007 

3007  K=K+1 

GO  TO  3003 

3008  IBA(I)»1 
IKILL( I )  «  J 
GO  TO  600 

f 3009  RN=RNG(NX) 

S1=T+(ASX( J)-ASM( J) )*RN+ASM(J) 
BME(ItJ)«Sl 
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IS1=5 

IS2=I 

IS3=J 

ASSIGN  3010  TO  NEXT 

GO  TO  700 

3010  IS1=2 
IS2  =  I 

ASSIGN  3011  TO  NEXT 
GO  TO  700 

3011  IS1=4 
IS2  =  I 

ASSIGN  600  TO  NEXT 
GO  TO  700 

3013  RN=RNG(NX) 

_S1=T+ ( ASX ( J ) -ASM ( J ) ) *RN+ASM ( J ) 
BME(I»J)«S1 
IS1*2 
IS2  =  I 
IS3=J 

ASSIGN  3014 
GO  TO  700 

3014  IS1=4 
ASSIGN  3008 
GO  TO  700 


TO  NEXT 


TO  NEXT 


C 
C 
C 


C 
C 

c 


ENGAGEMENT 

4000  I.  =  IE2 
J  =  IE3 
BME(I»J)=0. 
GO  TO  600 

FIRE  EVENT 


5000 


5001 
5002 
5003 
5004 
5005 
5006 
5020 

5007 
5008 

5010 
5011 


5012 


1  =  1 

J*I 

IF( 

IF( 

IF( 

IF 

IF( 

K=l 

IF( 

Sl  = 

GO 

IF( 

K  =  K 

GO 

Sl  = 

IS1 

IS2 

IS3 

ASS 

GO 

Sl  = 

GO 


E2 

E3 

T-BTPC(I»JtN) ) 5001* 600 » 600 

IBA(D)  600»  5002,  600 

MAL(J))  5003»  5010»  5003 

(MTRF(J))  5004.  5012,  5004 

IRUN)  5005»5013»5005 


BME( 
BME( 
TO  5 
K-NM 
+  1 

TO  5 
T+ST 
=  5 
*I 
=  J 
IGN 
TO  7 
T+ST 
TO  5 


UK))  5020»5007,5020 

I»K) 

Oil 

ISAR)  5008»5013»5008 

006 
EP 


600  TO  NEXT 

00 

EP 

Oil 
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5013  XA=BX( I»N)+BVX( I ,N )* ( T-BT ( I ,N) ) 
YA=BY( I »N)+BVY( I ,N)*(T-BT( I.N)  ) 
AfBS(I.N)**2  -  AVS(J)*»2 

B=  2.0  *  (BVX( I ,N)*(XA-AX(J) )+BVY< I .N) * ( YA-AY ( J ) ) ) 
C=  (XA-AX( J) )**2  +  (YA-AYU)  )**2+  ALTU.N)**2 
IF  (A)  5114,  5123*  5114 

5114  Z=  B**2-4.0*A*C 

IF  (Z)  600.  5115t  5115 

5115  TF1=  (-B  +  SQRTF(Z) )/(2.0*A) 
TF2«  (-B-SQRTF(Z) )/(2.0*A) 
IF  (TF1)  5116,600.5116 

5116  IF  (TF2)  5117.600.5117 

5117  IF  (TF1)  5120.5118,5118 

5118  IF  (TF2)  5119.5122.5122 

5119  TF=TF1 

GO  TO  5015 

5120  IF  (TF2)  600.  5121.  5121. 

5121  TF*TF2 

GO  TO  5015 

5122  IF  (TF1-TF2)  5119.  5121.  5121 

5123  TF  =  -C/B 
IF  (TF)600. 600. 5015 

5015  R=>  (TF  *  AVS(J))**2  -  ALT(I.N)**2 
R=  SQRTF(R) 

IF  HR-l.)-RMAX(J) )  5016.5016.5019 

5016  S1=T+TF 
IS1-3 
IS2*I 
IS3»J 

ASSIGN  5017  TO  NEXT 
GO  TO  700 

5017  MAL(J)=MAL( J)-l 
MTRFU)«MTRF(J)-1 
IBMEC( I»J)»IBMEC( I»J)+1 
BME(I.J)*S1 
RN»RNG(NX) 

S1=T+(RTX(J)-RTM( J))*RN+RTM(J) 
IS1«1 
IS2  =  I 
IS3»J 

5018  ASSIGN  600  TO  NEXT 
GO  TO  700 

5019  S1=T+STEP 
IS1-5 
IS2«I 
IS3»J 
GO  TO  5018 

C 

C      OUTPUT 

C 

'6000  IF(  OUTPUT)  9999.  6001.  6015 

6001  DO  6002  I=l.NBOMB 
ITOT(I)  ■  0 
DO  6002  J=1»NMISAR 

6002  IT0T(I)=  ITOT(I)  +  IBMEC(I.J) 
IF  (NFLAG-1)  6003.6004.6006 
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6003  IF(IRUN)  9999,6004,6006 

6004  WRITE  OUTPUT  TAPE  6,6005, N 

6005  FORMAT(/////»35H  UNCOORDINATED  ATTACK  RESULTS    REPtI3»//) 
GO  TO  6008 

6006  WRITE  OUTPUT  TAPE  6, 6007, N 

6007  F0RMAT(/////»33H  COORDINATED  ATTACK  RESULTS    REP»I3»//) 

6008  WRITE  OUTPUT  TAPE  6,  6009 

6009  FORMAT  (  45H  BMR  SMA1  SMA2  SMA3  TOTAL  A/K  XCOORD  YCOORD 

1   51HBRH   ALTITUDE  SPEED   TIME   TIME(l)  TIME(2)  TIME(3)  //) 
DO  6010  I=1,NB0MB 

6010  WRITE  OUTPUT  TAPE  6,6011,1,  ( I BMEC( I , J ) ♦ J=l »3 ) » ITOT ( I ) » 

1  IKILL( I ),BX(I,N) »BY( I,N) »BRH(I,N),ALT( I,N),BS( I,N) , 

2  BT( I»N), (BTMA( I,J,N) , J=l ,NMISAR ) 

6011  FORMAT  ( 14,  315.  16,  .14,  3F7.2.  F7.2,  F8.2,  F6.3,  F7.3.2F8.3) 
WRITE  OUTPUT  TAPE  6»  60115 

>0115  FORMAT  <///) 
60  TO  6015 

6015  I=N 

IF(  IRUN  )  9999,6016,6017 

6016  K=l 

GO    TO    6018 

6017  K=2 

6018  DO    6019    J=1,NMISAR 
CTOTSU»J»K)»0 
AM( I»J»K)*0 
DO    6019    L=1»NB0MB.    . 
B=IBMEC(L,J) 
CTOTS(I»J,K)    =CT0TS( I»J,K>    +B 

6019  CONTINUE 
DO    6021    L=1,NB0MB 
IFUKILL(D)    9999,6021,6020 

6020  J=IKILL(L) 
AM(I»J»K)=AM< I,J»K)+1.0 

6021  CONTINUE 
CT0T(I»K)  «.  0. 
AMT(I»K)  =0. 
DO  6022  J=1,NMISAR 
CTOT(I»K)=CTOT( I,K)+CTOTS( I,J,K) 

6022  AMT(I,K)=AMT( I,K)+AM< I,J»K) 
A=NMISAR 

ABAR(I,K)=AMT( I,K)/A 
CBAR(I,K)=CTOT( ItK)/A 
CVAR(I,K)  =0, 
AVAR(I,K)  *0» 
DO  6023  J=1,NMISAR 
CVAR(I,K)=CVAR( I ,K )+ ( CTOTS ( I » J, KJ-CBAR ( I ,K) )**2 

6023  AVARU»K)=AVAR( I»K)+(AM( I » J,K )-ABAR ( I >K ) )**2 
AVAR(I,K)=AVAR( I»K)/A 

CVAR(I,K)=CVAR( I,K)/A 
CSDU,K)»SQRTF(CVAR(  I»K)) 
ASD(I,K)=SQRTF(AVAR( I,K)) 
IF(N-NRPL)  7000,6024, 7000 
IF  (NFLAG-1)  60245,6025,6035 
IF  (IRUN)  6025,7000,6025 
K*l 
A«NRPL 
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DO    6027    J=1,NMI5AR 

DM( J»K)=0 

BM< J»K)=0 

DO    6026    I=1,NRPL 

DM(J»K)=DM( J>K)+CT0TS( I»JtK) 

6026  BM(J»K)=BM( J»K)+AM( I »J»K) 
DM( J»K)=DM( J»K)/A 

6027  BM(J»K)=BM( J»K)/A 
DMT(K)=0 
BMT(K)=0 

DO    6028    I=1,NRPL 

DMT ( K ) «DMT ( K ) +CTOT ( I  ♦  K  ) 

6028  BMT<K)=BMT(K>+AMTUtK) 
DMT(K)=DMT(K)/A 
BMT(K)«BMT(K)/A 

DO  6030  J=1»NMISAR 

DV(JtK)"0 

BV(JtK)  =  0 

DO  6029  I=1»NRPL 

DV(J»K)=DV( J»K)  +  (CTOTSU»J.K)-DM(J»K))**2 

6029  BV(J»K)=BV< J»K)+(AM( I » J »K)-BM{ J ,K ) )**2 
DV< J»K)=DV< J»K)/A 

6030  BV< J»K)=BV( J»K)/A 
DVT(K)«0 
BVT(K)»0 

DO  6031  I=1»NRPL 
DVT(K)=DVT(K)+(CTOT( I»K)-DMT(K) )**2 

6031  BVT(K)=BVT(K)+(AMT(I >K)-BMT(K) )  **2 
DVT(K)=DVT(K)/A 
BVT(K)=BVT(K)/A 
DO  6032  J=1,NMISAR 
DSD(J»K)=SQRTF(DV( J»K) ) 

6032  BSD(J»K)=SQRTF(BV(J»K) ) 
DSDT(K)*SQRTF(DVT(K) ) 
BSDT(K)*SQRTF(BVT(K) ) 
A=NBOMB 
PS(K)  =0. 
DO  6034  I=ltNRPL 
IF(AMT( I»K)-A)6034»6033»6034 

6033  PS(K)=PS(K)+1.0 

6034  CONTINUE 
A=NRPL 

PS(K)=PS(K)/A 
IF  (NFLAG-1)  60345»6036.6049 

60345  IF  (K-2)  6035 t 6036 »6035 

6035  K=2 
60  TO  60255 

6036  K=l 

6037  WRITE  OUTPUT  TAPE  6.6038 
,6038  F0RMAT( 31H1SUMMARY  UNCOORDINATED  RESULTS  //» 

1         19X»5HKILLS»51X»11HSH0TS  FIRED  /) 
GO  TO  6041 
6039  IF  (NFLAG-1)  60395»  60395»  60390 
60390  WRITE  OUTPUT  TAPE  6»  60405 
60405  FORMAT  (1H1) 
60395  WRITE  OUTPUT  TAPE  6t  6040 
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6040  FORMAT  ( 29H  SUMMARY  COORDINATED  RESULTS  //» 
1         19X,5HKILLS,51X,11HSH0TS  FIRED  /) 

6041  WRITE  OUTPUT  TAPE  6,6042 

6042  FORMAT  ( 34H  NBOMB  REP  MAI  MA2  MA3  TOTAL   MEAN,4X» 

1  9HVAR     SD,7X,20HMAl   MA2   MA3   TOTAL, 

2  19H   MEAN     VAR     SD/) 
WRITE  OUTPUT  TAPE  6,6043,  ( NBOMB, I , ( AM ( I ♦ J,K ) , J=l ,3 ) , 

1  AMT( I,K),ABAR( I  ,K) ,AVAR( I,K) ,ASD( I»K)  , 

2  (CTOTS( I,J,K) ,J=1,3) ,CTOT( I,K) ,CBAR( I,K) , 

3  CVAR( I,K) ,CSD( I,K) ,I=1,NRPL) 

6043  FORMAT  ( I  6, I4,3F4.0, F6«0,3F7«2 ,F8.0 ,2F5.0,F7 .0 »3F7.2 ) 
WRITE  OUTPUT  TAPE  6,  6044 

6044  F0RMAT(//////,9X,26HMA1     MA2     MA3     TOTAL»28X, 
1        26HMA1     MA2     MA3     TOTAL/) 

WRITE  OUTPUT  TAPE  6,6045,  (BM(J,K),  J=l,3)»   BMT ( K ) • 
1         (DM(J,K) tJ=l,3),DMT(K) 

6045  FORMAT(6H  MEAN  3F7.2  ,F8 .2 , 24X,3F7.2 »F8.2 ) 
WRITE  OUTPUT  TAPE  6,6046,  <BV(J,K)»  J=l,3)»  BVT ( K) » 

1         (DV( J»K),J=1»3) »DVT(K) 

6046  FORMAT  (6H  VAR   3F7. 2 , F8.2 ,24X ,3F7.2 ,F8 .2 ) 
WRITE  OUTPUT  TAPE  6,6047,  ( BSD( J ,K ) , J-l ,3 ) ♦  BSDT(K), 

1         (DSD( J,K) ,J=1.3) ,  DSDT(K) 

6047  FORMAT(6H  SD    3F7.2  ,F8.2 ,24X, 3F7.2 , F8. 2) 
WRITE  OUTPUT  TAPE  6,  6048,PS(K) 

6048  FORMAT  (///»  26H  PROBABLITY  OF  SURVIVAL  *  F7.2»/»1H1) 
IF  (NFLAG-1)  60485,7000,7000 

60485  IF  (K-2)  6049,7000,9999 

6049  K=2 
GO  TO  6039 

C 

C      REPEAT  AND  END 

C 

7000  NOL=0 
IF  (NFLAG-1)  900,7030,7030 

900  IF  (IRUN)  9999,7001,7030 

7001  IRUN=1 
ASSIGN  7020  TO  LAST 

7005  1=1 

7006  J=l 

7007  IF(BTMA(I,J,N) ) 7008, 7009,7008 

7008  S1=BTMA(I,J,N) 

151  =  5 

152  =  I 

153  =  J 

ASSIGN  7009  TO  NEXT 
GO  TO  700 

7009  IF(J-NMISAR)  7010,7011,9999 

7010  J=J+1 
GO  TO  7007 

7011  IF< I-NBOMB)7012, 7013, 9999 
'7012  I«I+1 

GO  TO  7006 
7013  GO  TO  LAST 
7020  DO  7021  1=1, NBOMB 

IBA(I)=0 

IKILL( I)»0 
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ITOT(I)=0 

DO  7021  J=1.NMISAR 
IBMECt I »J)=0 
BME( I»J)=0 

7021  CONTINUE 

DO  7022  J=1,NMISAR 
MAL(J)=MALO(J) 

7022  MTRF(J)=MTRFO(J) 
RN«RNG(0) 

RN=  RNG(NRl) 

IF  (NFLAG-1)  70225»70225»7060 
70225  IF  (IRUN)  9999»7060» 600 

7030  IF(N-NRPL)7031»7050»9999 

7031  IF(NREP)  9999.7032»7040 

7032  IF  (NFLAG-1)  70325 »70325»810 
810  IRUN=1 

GO  TO  70326 

70325  IRUN=0 

70326  NR1=NR1+NRDEL1 
NR2=NR2+NRDEL2 
N  =  N+1 

GO  TO  7020 

7040  IF  (NFLAG-1)  70405 »70405»820 
820  IRUN=1 

GO  TO  70406 

70405  IRUN=0 

70406  NR1<=NR1+NRDEL1 
N=N+1 

NBOMB=NBOMB-NBDEL 
ASSIGN  7041  TO  LAST 
GO  TO  7005 

7041  NR2*=NBOMB*4-1  +  IFE(N)+NR20  +  NRDEL2*(N-1) 
GO  TO  7053 

705  0  IF (NBOMB-NBMAX) 7051 »9888* 9888 

7051  NSOMB=NBOMB+NBDEL 

IF (NBOMB-NBMAX) 7052 f 7052*9888 

7052  NBOMB=NBOMB-NBDEL 

IF  (NFLAG-1)  70525»70525»830 
830  IRUN*1 

GO  TO  70526 

70525  IRUN'O 

70526  N=l 
NREP=1 

NR2=NBOMB*4-1+IFE(N)+NR20 
NR1=NR10 

ASSIGN  7053  TO  LAST 
GO  TO  7005 

7053  NBOMB*NBOMB+NBDEL 
GO  TO  7020 

7060  IF(NREP)9999»31»7061 

7061  I=NBOMB-NBDEL+l 
RN=  RNG(O) 

RN=  RNG(NR2) 

GO  TO  33 
9999  PRINT  9997 
9997  FORMAT  (  6H  ERROR) 

92 


9888  IF(  IEND)9994»9994»1 
9994  CONTINUE 

END 

FUNCTION  RNG(NR) 

CON  (K«30517578125» 

1  Kl=30517578125f 

2  K2*2000000000000000B» 

3  K3=7777400000000000B) 
IF  (NR)  10»10t20 

10  Kl  ■  K 

GO  TO  40 
20  DO  30  I«1»NR 

30  LDA(K) »MUI(K1)»SCL(K3) •STA(Kl) 
40  ALS(l) »ADD(K2).FAD(K2) 

RETURN 

END 

END 
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APPENDIX  IV 
Using  The  Program  at  The  U.  S.  Naval  Postgraduate  School 

For  the  FORTRAN -60  Compiler  and  the  CDC-1604  now  in  use  at 
the  U.  S.  Naval  Postgraduate  School  the  program  is  included  on  the 
library  tape.   The  program's  name  is  REGAME  and  using  the  program 
requires  only  that  the  user  submit  a  normal  FORTRAN-60  job 
consisting  of  a  simple  FORTRAN-60  program,  by  which  the  program 
REGAME  is  loaded  into  memory  by  the  statement  CALL  REGAME,  and  the 
necessary  sets  of  4  input  cards  as  described  in  the  Input  Section. 

Except  for  the  values  to  be  placed  on  the  input  cards  the  next 
page  of  this  appendix  is  an  example  of  the  FORTRAN-60  program 
needed  for  using  the  REGAME  program. 
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